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pThe  first  part  of  this  report  details  a  review  of  the  disperaion  of 
pollutants  by  rainfall,  with  particular  reference  to  inertial  and 
diffuaive  removal  by  precipitation,  inertial,  condenaational  and 
diffualva  removal  by  cloud  dropa,  and  the  influence  of  electrical 
forces  on  scavenging. 

The  second  half  of  the  report  describee  a  preliminary  study  of  the 
influanca  of  meteorological  parametera  on  the  size  distribution  of 
natural  aaroaol  particles,  with  particular  referencs  to  two  data  — 
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1 .INTRODUCTION 


Much  of  the  uork  described  in  this  report  constitutes  an 
extension  of  the  preliminary  studies  carried  out  in  the  three 
years  prior  to  December  1976  and  described  in  the  Final 
Technical  Report  under  contract  no.  DA-ERQ-124-74-G0D23. 

Two  main  laboratory  studies  undertaken  during  the  period  April 
1977  -  April  1980,  to  examine  the  scavenging  of  submicron 
particles  by  electrified  drops  and  the  influence  of  mixing  on 
the  evolution  of  droplet  spectra,  have  been  reported  in  Interim 
reports  and  uie  include  in  this  report  some  preliminary  studies 
of  the  dispersion  of  pollutants  by  rainfall. 

From  studies  made  in  the  UMI5T  laboratories  and  at  the 
field  station  of  the  Atmospheric  Physics  Research  Group  at 
Great  Dun  Fell,  Cumbria,  ue  describe  the  results  of  an  invest¬ 
igation  of  the  influence  of  meteorological  parameters  on  the 
size  distribution  of  natural  aerosol  particles  in  both  urban 
and  rural  environments. 
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2.  DISPERSION  OF  POLLUTANTS  BY  RAINFALL 

The  Earth  la  40%  covered  by  clouds:  this  figure  being 
arrived  at  by  average  albedo  measurements.  The  efficiency 
of  cloud  removal  processes  is  high  over  a  wide  range  of 
aerosol  particle  sizes  in  contrast  to  the  processes  occurring 
belou  cloud  uhich  mainly  occur  for  giant  particlea  (larger 
than  2  microns).  Uet  removal  is  usually  divided  into  rainout 
(in-cloud  processes)  and  uashout  (belou-cloud  processes). 
Evaporation  belou  cloud  base  uill  increase  the  concentrations 
of  the  constituents  and  the  presence  of  high  uind  uill  affrct 
the  number  of  drops  removed. 

2.1  Inertial  removal  by  precipitation 

Belou  cloudbase  sueepout  of  giant  aerosol  particles  in 
the  range  (2  to  10)  microns  is  the  major  removal  process. 
Although  this  process  may  occur  in  clouds,  the  majority  of  the 
time  the  droplet  is  in  flight  it  is  belou  cloud  base.  At  the 
5QQ  millibar  level,  raindrops  travel  at  50%  higher  velocities 
than  at  ground  level  (Beard  1976). 

Sueepout  involved  the  collection  of  particles  in  a  rain¬ 
drop's  path.  This  collection  depends  on  the  inertia  of  the 
particle  and  the  Reynolds  number  of  the  flou  around  the  drop. 

The  fraction  of  particles  captured  by  the  drop  to  the  particles 
beneath  the  drop  is  the  collection  efficiency  E. 

Assuming  unity  collection  efficiency,  the  area  suept  out 
depends  on  the  area  of  the  drop  (r  ),  uhilst  the  volume  of  the 
drop  (r3)  determines  the  net  concentration.  Thus  collection 
of  particles  belou  cloud  base  uould  result  in  K  «  —  uhere  K 
is  concentration  of  pollutant  in  drop  and  r  ia  radius  of  drop 
Evaporation  increases  the  relative  concentrations  in  the  smaller 
drops.  Turner  (1955)  found  for  freezing  clouds  that  K  ■  r”x 
uhere  r  is  radius  of  drop  and  K  the  concentration,  x  uas  found 

to  be  variable  and  he  attributed  this  result  to  evaporation  and 


collection  of  salt  oil ow  cloudbsse. 

According  to  Jungs  (1963)  the  net  concentration  in  the 

drops  (K")  ia  given  by 

K"  -  fK  ♦  K* 

H*  -  ~  Eb3N  (a)  da 
r 

f  is  a  term  to  represent  evaporation 
K'  is  the  washout  term 
a  ia  the  radius  of  an  aerosol  particle 
r  ia  the  radius  of  a  typical  raindrop 
E  la  the  collection  efficiency 
H  is  the  height  of  fall  below  cloudbaae 
K  represents  rainout 

N  ia  the  number  density  of  particles  radius  a. 

The  collection  efficiency  E  has  been  calculated  empirically 
by  many  workers  from  Langmuir  (1945)  onward,  e.g.  Beard  and 
Grover  (1974).  These  calculations  were  based  on  inertial 
deposition  on  a  drop  assuming  unity  coalescence  efficiency. 
However,  this  is  only  the  case  for  hygroscopic  particles  and 
some  of  these  will  already  have  been  removed  preferentially 
incloud  and  by  other  drops. 

Ranz  and  blong  (1952)  measured  the  collection  efficiencies 
for  spheres  as  a  function  of  Stokes'  number.  They  found  the 


simple  relation 


K  ■  Stk  *  *  (£> 


E  is  collection  efficiency 

Stk  ia  Stokes'  number 

p  ia  density  of  aerosol  particle 

V  ia  terminal  velocity  of  raindrop 

n  ia  viscosity  of  air 

Other  parameters  as  previously  defined. 

The  function  ~  does  not  change  very  much  with  radius.  However, 


variation  of  this  uith  aerosol  particle  radius  will  ba  a  step 
function.  Thus  particles  smaller  than  1.0  micron  will  hardly 
be  removed  but  particles  larger  than  thla  will  all  be  removed. 
Storebo  and  Dingle  (1974)  looked  at  pollution  removal  theoretically 
and  came  to  the  conclusion  that  hygroscopic  particles  were  easily 
removed.  Increased  windepeed  wee  shown  to  give  a  decade's  differ- 
ence  in  concentrations  at  25  kilometres,  probably  because  of 
increased  condensation  rate  and  leas  removal  of  drops. 


2.2  Removal  by  diffusion  to  precipitation 

Diffusive  removal  by  precipitation  1b  treated  by  tha  advection 
diffusion  equation 


dn 

dr 


on  ivc  . .  _ 

— -  «  u  n  -  v  n 


n  ia  number  of  particles  per  unit  volume 
v  la  velocity  of  fluid 
D  is  diffusion  constant 

(See  Twomey  1977).  The  collection  efficiency  1b  specified  by 


the  Reynolds  number 


prV 


and  peclet  number 


2rV 


p  la  Penalty  of  fluid 
r  la  radius  of  raindrop 
V  is  terminal  velocity  of  drop 
D  la  diffusion  constant 

Frledlardsr  derived  the  expression  (1967) 


E  *  4  Pe”1  •  (1  +  0.2  Re1/6  Pe1/3) 

E  is  diffusion  collection  efficiency 
Pe  la  Peclet  number 
Re  is  Reynolds  numbar 

The  expression  (0.4  Re1^  Pe1^3)  varies  weakly  with  Reynolds 


number  and  ia  negligibly  small,  hence  it  can  be  deduced  that 


bo  that  for  an  order  of  magnitude  change  in  drop  radius  (e.g. 
from  drizzle  to  raindropa)  the  efficiency  is  changed  by  approxi¬ 
mately  a  factor  of  IQ.  Although  the  efficiencies  of  removal 
appear  very  low,  the  numbers  of  particles  in  the  lower  size 
ranges  means  that  considerable  amounts  of  aerosol  will  be 
removed,  e.g.  below  Q.1  microns.  For  raindrops  and  particles 
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□.4  microns  the  efficiency  is  IQ  and  this  is  very  close  to  the 
experimental  limits  of  measurement. 

The  intermediate  range  between  diffusion  and  inertial 
removal  (Q.1  to  1,Q  microns)  has  been  investigated  by  Kerker 
and  Hampl  (1974)  who  found  experimentally  that  Zebel'a  expression 
was  valid: 

E  »  1.68  Pe”2/^  which  is  in  the  range  1Q_if  to  10”^. 

2.3  Inertial  removal  in  cloud 

The  terminal  velocity  of  cloud  drops  is  small  (1cm  per 
second).  Falling  cloud  drops  are  net  very  effective  at  scavenging 
smaller  particles  even  though  coalescence  into  raindrops  involves 
a  very  large  number  (about  1  million)  of  cloud  drops. 

In  the  case  of  clouds  which  comprise  ice  crystals,  e.g. 
cumuliform  types,  the  ice  crystals  may  act  as  scavengers  of 
aerosol  particles.  Pitter  (1977)  has  numerically  modelled  the 
scavenging  of  aerosol  particles  by  thin  plates  under  atmospheric 
conditions  of  -18°C,  400  millibars,  ice  crystal  lengths  being 
103  -  366  microns.  The  model  included  hydrodynamic,  gravitational 
and  electrostatic  forces.  He  concluded  that  for  103  microns 
crystals  ,  scavenging  efficiencies  of  0.05  to  0.4  for  aerosol 
particles  from  0.5  to  10  microns.  He  found  that  for  uncharged  ice 
crystals,  there  was  a  cut  off  at  5  microns.  This  also  applied 
to  replusive  charges.  Sood  and  Jackoon  (1970)  investigated 
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experimentally  snow  scavenging  far  particle  sizes  0.25  microns, 

□.5  microns  and  came  to  the  conclusion  that  there  is  a  pronounced 
minimum  in  scavenging  efficiency  at  about  0.4  microns*  Scavenging 
efficiencies  of  the  order  0.6  to  2%  were  found  for  sodium  chloride 
aerosol  particle  sizes  0.1  to  0.16  microns. 
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2.4  Removal  by  nucleation 

Cloud  drops  will  begin  to  grout  when  a  critical  super- 
saturation  has  been  maintained.  A  typical  hygroscopic  particle 
acting  as  condensation  nucleus  mill  undergo  considerable  dilution 
as  the  drop  grouts  from  critical  radius  (Q.02  microns)  to  a 
typical  cloud  drop  size  of  10  microns.  Katz  and  Kocmond  (1973) 
concluded  that  cloud  condensation  nuclei  may  be  considerably 
larger.  According  to  Rogers  (1976),  typical  droplet  number  den¬ 
sities  lie  betuieen  50  and  200  per  cc  although  aerosol  size  dist¬ 
ributions  are  very  variable  only  a  small  number  of  particles 
serve  as  cloud  condensation  nuclei,  e.g.  for  the  Continental  size 
distribution  of  Junge  (1963)  10%  serve  as  nuclei.  Following 
nucleation,  the  resulting  cloud  drop  has  sufficient  inertia  to  grow 
by  condensation  and  accretion  to  be  removed  gravitationally.  If 
a  cloud  formed  by  condensation  of  the  largest  drops  on  the  largest 
cloud  nuclei  the  subsequent  growth  by  coalescence  to  raindrop 
size  would  occur  relatively  quickly.  Turner  (1955)  applied  this 
argument  to  non-freezing  clouds  in  Hawaii  to  account  for  the 
minima  found  in  his  results.  Coalescence  is  still  an  important 
growth  mechanism  for  cloud  drops  in  warm  clouds  with  freezing 
tops  (see  Rogers,  1976).  Takahashi  (1976)  developed  a  computer 
model  for  warm  clouds  which  correlated  with  Turner's  results  for 
Orographic  Hawaiian  rain. 

2.5  Diffusion  to  cloud  drops 

Brpwnltn  diffusion  of  primarily  Altktn  par t i else  to  cloud  ..jm 


T 


less  than  0.05  microns.  Junge  (1963)  has  shown  from  the  Smoluchowsky 
formula  that: 


cl?  “  <*ndrcnc  9ive8  8  half  life  of 

T  *  0.69/4n  d  r_n_  where 
c  c 

n  *  concentration  of  Aitken  particles 
d  *  diffusion  constant  of  Aitken  particles 
rQ  ■  radius  of  cloud  drops 
nQ  «  number  of  cloud  drops 
t  a  time 

and  he  goes  on  to  show  for  r_  a  10  microns  and  n„  a  200  per  cc 

c  c 

(bJeickmann  and  Kampe  1953)  corresponding  to  a  liquid  water  content 
of  0.84gm  per  cubic  metre 
r  microns  0.02  0.03  0.1 

T  hours  0.64  3.8  38 


The  estimate  will  not  change  much  if  charges  are  considered. 
This  treatment  is  justified  since  the  separation  of  drops  is 
about  1mm  ao  that  there  is  no  interaction  of  the  drops.  The 
raindrops  farmed  by  coalescence  of  these  cloud  droplets  will 
all  have  the  same  concentrations  of  contaminants.  Goldsmith, 
Delafield  and  Cox  (1963)  have  shown  that  this  process  is  very 
efficient  for  particles  smaller  than  0.02  microns. 


3°ptpc 


where 


lilgQ  is  the  washout  coefficent  defined  by 

w  specific  activity  of  rain 

UBD  *  specific  activity  of  air 

in  the  case  of  the  radioactive  scavenging  studies.  Electric 

charges  were  assumed  to  have  negligible  effect  in  the  presence 

of  ionising  radiation. 

0  is  the  diffusion  constant 


p  is  the  density  of  aerosol 
Pa  ia  the  density  of  eir 


p  is  the  density  of  uater 
a  is  the  radius  of  aerosol  particles 
t  is  time 

The  time  scale  ia  quite  aignificant  e.g.  for  a  duration  of 
IQ  hours  all  the  particles  radius  Q.Q2  microns  and  most  of 
□•□3  microns  will  be  removed.  The  flux  of  aerosol  particles 
along  the  vapour  gradient  (Facy  effect)  may  be  important  for 
particles  smaller  than  0.1  microns.  Goldsmith  (1961)  obtained 
lou  rates  of  attachment,  though  this  took  no  account  of  electric 
charges. 

Phoretic  forces,  produced  by  ordered  fluxes  superimposed 
on  the  relatively  disordered  motion  in  a  gas  are  thought  to 
contribute  very  little  to  scavenging.  Dif f uaiophoreais  is 
caused  by  uater  vapour  produced  in  connection  uith  evaporation 
and  condensation.  Thermophoresia  occurs  similarly.  Houever,  the 
integrated  flux  ia  negligibly  small  for  the  condensation 
processes  only  last  for  10  seconds  of  a  typical  cloud  drop's 
life.  For  particles  in  the  sub-micron  range,  the  effect  uas 
shaun  to  be  negligibly  efficient  by  Goldmaith,  Delafield  and 
Cox  (1963),  typically  0.1  -  0.5J4.  Slinn  and  Hales  (1970)  have 
investigated  the  role  of  phoretic  forces  in  the  scavenging  of 
particulates  in  the  size  range  0.01  to  1.0  microns  and  found 
efficiencies  in  the  range  Q.00Q5  to  0.01  for  thermophoresis 
for  raindrops  0.1mm  to  1mm.  Dif f usiophoresis  uas  about  one 
third  as  efficient.  These  processes  may  occur  uhile  droplets 
are  grouing  before  coalescence  becomes  the  main  removal  process. 
Using  radioactive  tracers,  Perkins,  Thomas,  Young  and  Scott  (1970) 
found  that  light  rain  takes  longer  to  develop  than  heavy  rain. 

They  also  found  that  rainout  uas  initially  high,  then  decreased 
during  a  thunderstorm  and  touards  the  end  of  the  thunderstorm 
uould  increase  again.  They  concluded  that  crops  may  be  subject 
to  several  cycles  and  that  in  heavy  rain,  dif f usiophoresis  may 


be  significant.  Greenfield  (1957)  has  investigated  radioactive 
fall-out  using  a  model.  It  mas  found  that  there  uas  a  pronounced 
minimum  of  particles  scavenged  in  the  range  0.1  to  1.0  microns 
but  belou  0.01  microns  snd  above  10  microns  nearly  all  the 
particles  ue re  removed.  For  particles  2  microns  radius  rainfall 
rate  2.5mm  per  hour,  90%  of  the  total  activity  uas  removed,  the 
cloud  drops  in  each  case  uere  20  microns.  Particles  belou  1  micron 
mere  assumed  to  mix  uith  the  uater  cloud  before  rain  started. 

These  could  then  be  scavenged  by  cloud  droplets  by  a  coagulation 
process  and  hence  removed  by  the  rain. 

2.6  The  effect  of  electric  forces  in  scavenging 

Dang,  Grover  and  Pruppacher  (1976)  have  modelled  the 
collision  efficiencies  cf  raindrops  and  aerosol  particles  for 
different  radii,  taking  electric  charges  and  different  humidities 
into  account.  Drop  sizes  6.6,  18.6,  30,  UZ,  72,  106,  173  microns 
uere  used.  The  electric  charges  uere  claimed  to  be  realistic. 

The  efficiencies  decreased  quite  markedly  from  about  one  to 
one  five-hundredth  for  particle  radius  0.1  microns.  The  effect 
of  electric  charge  uas  quite  marked,  the  tendency  being  to 
increase  collection  efficiency  up  to  tuo  orders  of  magnitude  over 
the  range  0.005  to  1.0  microns,  the  resultant  efficiencies  being 
as  high  as  one  for  8.6  microns  drops. but  about  0.005  for  173 
micron  drops.  The  effect  of  electric  forces  on  scavenging  has 
been  investigated  by  Dayan  and  Gallitz  (1975)  experimentally  and 
they  found  that  the  collection  efficiency  uas  increased  by  1 
to  2  orders  of  magnitude  over  uncharged  drops.  May  (1958)  in¬ 
vestigated  the  uashout  of  lycopodium  spores  by  rain  and  found 
charges  of  about  6  *  10-i*esu.  Smith  (1955)  shoued  that  charging 
on  drops  varies  from  10  to  5  *  10  esu:  Adam  and  Semonin  (1970) 
looked  at  submicron  scavenging  and  found  that  for  a  charged  drop 
radius  200  microns,  collection  efficiency  uas  about  10%,  though 
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these  cherges  may  be  high.  For  larger  drops  end  uncharged  drops 
this  use  negligible.  The  experimental  work  mas  carried  out 
using  bacillus  spores.  Charged  ice  crystals  ere  also  thought 
to  contribute  to  submicron  scavenging. 


3.  THE  INFLUENCE  OF  METEOROLOGICAL  PARAMETERS  ON  THE  51 


DISTRIBUTION  OF  NATURAL  AEROSOL  PARTICLES 


3.1  Introduction 

Ground-based  measurements  of  the  concentration  and  size 
distribution  of  aerosol  particles  in  the  range  0.5  to  10|jm 
diameter  were  made  at  both  Manchester  and  Great  Dun  Fell,  using 
a  Royco  optical  particle  counter.  They  are  described  and  in¬ 
terpreted  in  terms  of  the  prevailing  meteorological  conditions. 

A  considerable  amount  of  data  exists  defining  the  aerosol 
particle  spectrum  for  a  wide  range  of  localities  and  meteorological 
conditions,  though  the  manner  in  which  the  particulate  spectrum 
has  been  influenced  by  the  synoptic  meteorological  situation 
has  usually  received  little  attention.  However,  some  studies 
have  investigated  the  relationship  between  previling  meteorological 
conditions  and  bulk  concentrations  af  pollution. 

In  an  attempt  to  encompass  both  of  these  approaches  we  have 
undertaken  preliminary  studies  of  the  influence  of  meteorological 
parameters  on  the  size  distribution  of  large/glant  aerosol 
particles  (O-Sym  <  particle  diameter  <  lO^m)  for  prolonged 
sampling  periods.  Measurements  have  been  made  at  two  widely 
different  locations,  one  urban  with  ita  concomitant  anthropogenic 
sources  of' pollution,  and  one  remote,  rural,  high  level.  The 
meteorological  data  used  was  routine  measurements  available  from 
the  national  weather  service. 

3.2  Measurement  Sites.  Experimental  Estimates  and  Data  Sources 


The  city  of  Manchester  is  situated  in  a  basin  surrounded  by 
the  Pennine  Hilla  to  the  north,  east  and  south  reaching  to 
~500m  aal  in  the  north  and  east.  To  the  west  the  land  la 
comparatively  flat  encompassing  the  Lancashire  plain  and  Mersey 
estuary.  Measurements  of  particle  size  distribution  were  made 
from  the  roof  of  the  main  building  of  UMIST  campus  which  ia 


adjacent  to  the  city  centre.  By  contract  the  UM1ST  field 
station  on  Great  Dun  Fell  (GDF)  in  Cumbria  la  847m  aal  and 
located  on  an  laolated  peak  on  the  Pennine  ridge  above  the 
Eden  Valley  which  runa  northueat/aoutheaat. 

The  particle  alze  distribution  in  the  alze  range  0.5** 

10^m  waa  monitored  uaing  a  Royco  model  225  optical  particle 
counter  with  model  507  digital  diaplay  and  autoacan  facility 
and  model  127  digital  printer.  .  The  autoecan  facility  permlta 
sequential  monitoring  of  particle  concentrations  in  the  five 
size  ranges  0.5-0. 7,  0.7-1. 4,  1.4-3,  3-5  and  5-10^m  diameter. 

Prior  to  its  initial  use  the  counter  was  calibrated  uaing 
polystyrene  latex  spheres. 

Meteorological  data  waa  obtained  from  the  Mancheater 
Weather  Centre  (MUC)  from  their  detailed  logs  together  uith 
upper  air  data  taken  from  the  Dally  Aerological  Record  (DAR) 
for  Aughton,  about  50km  west  of  Manchester  and  other  sonde 
atationa  around  Britain  when  appropriate.  By  uaing  the 
6-hourly  synoptic  charts  in  the  Daily  Weather  Report  (DWR) 
trajectories,  of  different  airinesses  could  be  obtained  as  well 
as  information  about  the  location  of  frontal  systems.  Because 
there  is  no  radiosonde  station  near  GDF  a  combination  of 
aacenta  at.  Shanwell,  Aughton  and  Long  Kesh  together  with  a 
consideration  of  air  trajectory,  can  be  used  to  infer  the 
structure  of  the  atmosphere  in  the  vicinity  of  the  mountain. 

Aerosol  data  for  the  periods  9-27  February  1978  in  Manchester 
and  10-16  September  1977  at  Great  Dun  Fell  are  diacuaaed  herein. 

The  period  10-21  February  1978  was  dominated  by  a  polar,  continental 
airatream  and  22-27  February  by  a  mainly  maritime  regime.  The 
GDF  data  ia  influenced  by  a  basically  unstable  polar  maritime 
flow  with  a  strong  subsidence  inversion  at  about  900mb  close 
to  the  level  of  the  Great  Dun  Fell  summit. 


3.3  Results  and  Discussion 


(i)  Period  10-16  February  1978  (Manchester  data) 

'  Throughout  this  week  e  polar  continental  elretreero  pre¬ 
vailed  although  three  subtle  changes  of  airmaea  can  be 
detected.  From  00h  10/2  until  OOh  12/2  the  alratream  had 
Its  source  region  In  Northern  Finland  (airiness  1),  thereafter 

I 

the  source  was  mainly  Western  Russia  (airiness  2)  with  the 
airflow  arriving  at  Manchester  from  a  mainly  N-NE  direction. 

As  a  depression  moved  north  from  the  Mediterranean  into  the 
North  Bee  late  on  the  14th,  Manchester  came  under  its  cyclonic 
influence.  Air  circulating  around  this  system  came  mainly  from 
the  near  continent  and  arrived  from  a  mainly  east  to  south¬ 
easterly  direction  (airmaea  3).  The  trajectories  of  these 
airmasses  are  illustrated  in  Figure  1. 

From  the  diagrams  of  the  aerosol  counts  against  time 
(Figs. 2  &  3)  it  is  apparent  that  the  long  range  trajectory 
la  only  of  secondary  importance  in  determining  the  levels  of 
various  particle  sizes  on  these  occasions.  Shorter  range 
trajectories  appear  to  be  much  more  significant.  Broadly 
the  counts  in  0.35-0. 7pm  radius  range  show  an  increase  through¬ 
out  the  week  (10/2-16/2)  especially  after  13/2  end  this  is 
well  correlated  with  the  change  in  local  wind  direction  in 
Manchester  rather  than  changes  in  the  source  of  the  air. 

Early  in  the  week  the  direction  is  mainly  010°,  and  has  backed 
to  320°  on  the  13th.  Thereafter  there  is  a  steady  veer  through 
North  and  East  to  110°.  Figure  4  illustrates  the  changes  in 
count  of  the  0.35-0. 7^  range  with  wind  direction.  This 
particular  size  range  of  particle  has  its  highest  concentration 
when  the  air  has  come  from  the  industrial  West  Midlands  of 
England.  However  giant  particles,  eg,  2.5-5.0um  radius  range 
do  not  show  this  effect  and  this  would  suggest  e  local  source 


for  those  particular  aerosols,  reflecting  the  shorter  airborne 
lifetimes  of  the  larger  particles. 


The  stability  of  the  boundary  layer  proves  very  important 
in  determining  aeroaol  concentrations  by  affecting  the  depth 
to  which  turbulent  mixing  either  by  convection  or  shear  can 
occur.  Values  for  the  Richardson  number  have  been  computed 
from  radiosonde  data  at  Aughton  for  the  relevant  periods:  the 
Richardson  number  (Ri)  la  given  by 


where  g  ■  gravitational  acceleration 

o’  *  mean  potential  temperature  of  the  layer 

ae 

s  potential  temperature  gradient  of  the  layer 


«  shear  in  the  particular  lByer. 

This  parameter  effectively  measures  the  relative  production 
ratea  of  turbulent  kinetic  energy  by  buoyancy  (free  convection) 
to  ahear  (farced  convection).  Ri  >  0  ia  characteristic  of 
stably  stratified  layers  while  Ri  <  0  implies  instability  and 
thus  buoyant  motions  will  be  present. 

Figures  5  and  6  show  the  relationship  between  the  Richardson 
number  and  the  aeroaol  counts  for  both  3.5-0.?|,im  and  2.5+5.Qtjm 
size  ranges.  The  data  was  plotted  for  airmaBS  types  2  and  3. 

The  highest  counts  for  small  particles  are  found  to  coincide 
with  periods  of  strong  inversion  conditions  (large  and  positive 
Ri)  aa  emphasized  by  data  on  14/2  where  a  very  intense  nocturnal 
inversion  is  present.  During  this  particular  week  the  stability 
of  the  lowest ’ layers  decreased  markedly  during  the  day  (eg  on 
14th  aa  shown  by  the  Aughton  ascent  data  in  Figure  7)  aa  aolar 
heating  initiated  convection  and  subsequently  concentrations  of 
small  particles  decreased  sharply.  The  14th  ia  a  prime  example 


of  thlo  offset  whtn  following  diaptrool  of  tho  inversion 


eonvsetion  was  established  up  to  750mb  and  snow  showers  were 
reported.  For  the  large  particle  sizes  (2.5pm  radius)  there 
is  a  negative  correlation  between  concentration  and  Richardaon 
nuaber.  The  implication  la  that  during  periods  of  large  and 
positive  Richardaon  numbers  (high  stability  and  thus  little 
turbulence)  the  gravitational  settling  of  these  particles  ia 
important  and  sufficient  to  account  for  the  low  concentrations 
of  these  particles  in  the  early  mornings  of  these  days.  Further 
reference  to  the  diurnal  variation  of  the  particle  aizea  will 
ba  made  later  on  in  this  report. 

Cii)  Period  21-27  February  1978  (Manchester  Data) 

This  period  is  rather  complex  meteorologically  with  a 
fundamental  change  of  synoptic  type  early  in  the  period.  From 
OQh  21/2  to  Q6h  22/2  there  is  a  gradual  encroachment  of  Atlantic 
air  behind  a  aeriea  of  weak  occlusions  and  true  maritime  air 
does  not  arrive  until  the  24th.  Some  of  the  airatream  trajectories 
are  illustrated  in  Figure  8.  After  the  24th  a  aeriea  of  troughs 
and  fronts  in  the  generally  5-SW  flow  brought  periods  of  rain 
or  showers.  Generally  however  the  concentrations  of  particles 
in  all  size  ranges  sre  little  chenged  from  those  of  s  week 
earlier  in  the  continental  regime,  again  emphasising  the  import¬ 
ance  of  ehort/medium  range  trajectory. 

Aa  with  the  previoua  week' a  data,  the  concentration  of 
aeroaola  is  closely  aasociated  with  stability  for  a  given  air- 
mass.  An  intense  low-level  frontal  inversion  provides  an 
extremely  effective  barrier  for  the  dispersion  of  particulate 
matter  around  QOh  on  22/2  while  a  nocturnal  inversion  on  the 
27th  ia  aleo  responsible  for  high  concentrations  of  email 
particles.  Ae  for  the  14th,  the  27th  showe  similar  characteristlca 
with  daep  convactlon  occuring  ones  the  inversion  haa  dispersed 


eventually  producing  rain  ehowere.  Tha  axtraaaly  low  con- 
eentratlona  of  tha  24th  ara  aimilarly  aaaociated  with  a  very 
deap  convactlva  Mixing  layer,  and  thua  large  turbulence 
giving  riaa  to  coaparatlvely  large  concentratione  of  giant 
particles. 

Froa  OOh  on  23/2  the  analler  particle .concentration  rieea 
and  la  closely  aaaociated  with  a  veering  of  the  wind  froa  E  to  S, 
with  the  boundary  layer  stability  reaalning  easentially  un¬ 
changed,  thus  confiraing  the  elgniflcance  of  aediua  range  air 
trajectories;  the  industrial  Midlands  la  the  likely  eource  again. 

Figure  9  illustrates  the  diurnal  variation  of  the  particle 
concentrations  during  the  period  of  data  at  Manchester.  For 
particle  alzes  froa  0.25^a  up  to  1.5^0  radius  there  is  little 
evidence  of  any  diurnal  variation  but  there  are  strong  variations 
for  particles  upwards  of  2.0^m  radius.  The  reason  for  this  la 
the  diurnal  variation  of  turbulence  in  the  boundary  layer,  and 
therefore  Richardson  nuaber.  As  convection  becoaea  established 
in  the  early  morning  after  about  0700hrs,  the  turbulence  in  the 
boundary  layer  increases  (Ri  decreasing)  eventually  reaching 
a  maximum  intensity  at  about  150Qhrs.  Thereafter  convective 
turbulence  declines  and  the  ability  of  the  turbulent  fluctuations 
to  retain  giant  particles  aloft  is  diminished.  Because  these 
giant  particles  have  sufficiently  large  terminal  velocities 

_  A 

(~  1.0cm  a  ), gravitational  settling  becomes  an  Important 
process  in  determining  their  concentration  and  a  minimum  con¬ 
centration  ie  observed  around  OfiOOhra  which  ia  the  time  of 
minimum  turbulence,  le,  maximum  Richardson  number.  In  summery 
then,  the  diurnal  variation  of  stability  in  the  lowest  layers 
will  be  the  moat  fundamental  factor  in  determining  the  concent¬ 
ration  of  giant  particles. 

To  teat  this  conclusion,  correlation  coefficients  of 
Molt  ties  concentre tlone  against  varloue  aeteorologieel 


parameters  were  calculated  and  the  reaulta  are  preaented 
In  Table  1.  Illustrative  graphs  of  count  against  mateorologieal 
parameter  are  preaented  in  Figures  IQ  and  11.  Table  2  examines 
the  significance  of  these  calculations.  Strong  correlations 
between  wind  speed,  temperature  and  Ri  for  the  giant  particles 
provide  more  evidence  of  the  effect  of  the  diurnal  variation 
of  turbulence.  For  the  amall  particles  there  is  a  strong 
positive  correlation  with  the  Richardson  number,  but  the 
large  positive  values  of  Ri  are  generally  associated  with 
intense  low  level  inversions  either  frontal  or  as  la  more 
usual,  to  radiation  cooling.  Lower  values  of  Ri  are  associated 
with  greater  mixing  depths  in  convective  situations.  Elsom 
(1979)  provides  evidence  of  low  level  trapping  by  a  temperature 
inversion  which  because  of  its  longevity  gives  rise  to  anomalousl 
high  levels  of  pollution  in  Manchester. 


(iii)  Period  16-20  September  1977  (SDF  data) 

A  weak  trailing  cold  front,  associated  with  a  deep  depreaion 
over  Uestern  Russia,  lying  from  Essex  to  Northern  Ireland  at 
□Oh  16/9  moved  away  slowly  South-Ueat  and  had  cleared  the 
British  Isles  by  OOh  17/9.  Meanwhile  a  large  anticyclone  situated' 
over  Iceland  at  IZh  15/9  had  moved  slowly  South-east  to  lie 
between  the  Faeroe  Ielends  and  Shetlands  by  IZh  17/9  (Figure  12). 

h 

Thereafter  the  eynoptic  situation  remained  eeaentlally  unchanged 
with  the  anticyclone  becoming  stationary  and  keeping  much  the 
same  intensity.  The  effect  of  theee  synoptic  changes  was  that 

1 

the  winds  remained  north  -easterly  throughout  the  period  although  j 

t 

subtle  changes  of  trajectory  occurred.  The  alrstream  on  the 
16/9  and  17/9  had  its  origin  over  Greenland  and  had  reachad  60F 
around  the  eastern  and  south-eastern  flanks  of  the  anticyclone. 

By  the  end  of  the  period  as  the  anticyclone  had  slipped  south- 

_ _ 


eastwards,  tha  alrstresm  waa  coming  from  an  Atlantic  origin. 

Aa  theaa  air  maaaea  caaa  aouth  acroaa  the  North  Sea,  heating  by 
the  bceane  had  cauaed  them  to  become  unetable  in  their  loweat 
lay era  which  were  capped  by  a  atrong  aubaidenca  Inversion 
aeeociated  with  the  anticyclone*  The  trajectories  are  shown 
in  Figure  13* 

Obeervationa  of  a  full  atratocumulua  cloud  cover  at  the 
eaat  coaatal  atationa  of  Scotland  and  Northern  England  are 
characteriatic  of  thia  effect,  thua  particulate  matter  will  be 
well  mixed  in  the  leyer  up  to  the  lnvereion* 

A  graph  of  particle  concentration  against  time  (Figure  14) 
illustrates  large  variations  which  were  of  similar  magnitude 
and  were  in  phase  for  all  the  size  categories*  As  stated  pre¬ 
viously  the  meteorological  conditions  are  essentially  unchanged 
and  thua  changes  of  alrmaaa  cannot  account  for  the  large  fluct¬ 
uations, eg,  OOh  19/9  to  14h  19/9*  However,  changes  in  the 
height  of  the  eubeidence  inversion  are  very  noticeable  from 
radiosonde  aacenta  at  Shanwall  and  Aughton.  From  a  consideration 
of  air  trajectory  and  dlatanca  from  the  centre  of  the  anticyclone 
it  la  thought  that  theaa  two  stations  provide  a  representative 
sounding  of  the  air  in  tha  vicinity  of  Great  Dun  Fell.  Since 
the  meteorological  data  ascents  are  only  mode  at  OOh  and  12h 
during  the  day  so  the  eontopt  of  mixing  height  haa  to  be  employed 
to  extrapolate  valuea  of  the  inversion  height  at  times  between 
the  etendard  observations*  The  idee  of  mixing  height  was  proposed 
by  Holzworth  (1967)  and  the  maximum  mixing  height  ia  derived  as 
the  intersection  of  the  dry  ediabet  extanded  from  the  maximum 
temperature  to  the  aounding  at  12h.  The  minimum  mixing  height 
is  the  intersection  of  the  dry  adiebet  from  the  minimum  temper¬ 
ature  (plus  5*0°C  in  urban  areas  to  account  for  heat  island 
effects)  to  the  sounding  at  OOh. 


Figure 


height  with  time.  The  inference  to  be  drawn  ie  that  the 

V 

periods  of  very  low  particle  count  are  uh.en  the  invereion  falls 
below  the  mountain  top  and  the  eite  ie  in  a  region  where  the 
air  has  come  from  high  up  in  the  troposphere  by  subsidence. 

Higher  particle  concentrations  are  always  observed  during  the 
day  and  it  ie  likely  that  turbulent  mixing  caused  by  surface 
convection  from  the  velley  pushes  the  inversion  above  the 
mountain  top  with  the  inversion  providing  an  effective  bar 
against  further  dispersal.  Consequently  when  convection  dies 
away  again  in  the  evening  the  inversion  falls  below  the  mountain 
top  and  the  site  once  again  enters  the  subsidence  zone.  This 
effect  occurs  three  times  during  the  period.  (The  sampling  period 
wae  terminated  when  the  local  standby  diesel  generator  waa  activ- 

I 

ated,  resulting  in  the  sudden  increase  in  particle  concentrations  | 

i 

on  20/2)  Aacenta  at  both  Shanwell  and  Aughton  indicate  a  well- 
mixed  layer  below  the  inversion  during  the  day  with  a  auper- 
adiabatic  lapse  rate  established  near  the  ground  and 
a  dry  adiabatic  lapse  above  ity  up  to  the  inversion. 

Throughout  the  period  the  winds  at  9D0m  at  both  Shanwell 
and  Aughton  remain  less  than  20kts  from  the  North-east.  If  the 
wind  is  strong  it  is  observed  that  inversion  layers  can  be 
lifted  bodily  over  GDF  because  the  flow  has  sufficient  energy 
to  do  the  work  required.  For  lighter  airflows  this  ie  not 
necessarily  the  case  and  the  inversion  can  fall  below  mountain 

i 

top  level  particularly  when  a  strong  subsidence  inversion  if 
present.  Thus  it  may  be  deduced  that  the  variation  in  the  height 
of  the  Inversion  layer  ie  sufficient  to  account  for  the  fluct¬ 
uations  obaerved  in  the  aerosol  spectra. 


TABLE  1 

Correlation  Coefficients  for  Manchester  Data 


Partiela  Biza 
(li*  radius) 

0.25-0.35 


2. 5-5.8 


0.35-0.7 

2.5- 5.0 
0.25-0.35 

1. 5- 2.5 


Parameter 

Bind  apeed 
Relative  Humidity 
Temperature 

Wind  apeed 
Relative  Humidity 
Temperature 

Richardaon  number 

Richardson  number 

Rainfall  amount 

Rainfall  amount 


Correlation 
Coefficient  (r) 

♦  0.42 

-  0.27 
-i-  0.33 

♦  0.87 

-  0.88 

♦  0.95 

♦  0.94 
-0.85 
-0.17 
-0.43 


2. 5-5.0 


Rainfall  amount 


0.47 
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Airetream  trajectories  to  Manchester  in  the  period  10th- 
16th  February  1976. 

1  10.2.78;  2  13.2.78;  3  16.2.78 

Particle  concentration  N  cm*''*  versus  time  in  Manchester 


from  10th-16th  February  1978  for  channels  1-5. 

radius  range  covered  by  channel 
channel  1  0.25  -  0.35u» 


channel  2 
channel  3 
channel  4 
channel  5 


0.35  -  0.7  u® 

0.7  -  1.5  um 

1.5  “  2.5  ub 

2.5  -  5.0  pm 

Particle  concentration  N  cm"3  veraua  time  in  Manchester 
from  21at-27th  February  1978.  Size  ranges  covered  are  thoi 
for  Figure  2. 

Scatter  diagram  of  wind  direction  8  veraua  percentage 
deviation  from  the  average  concentration  D  for  channel  2. 
The  average  concentration  is  taken  for  the  period  10th- 
16th  February  1978. 

Richardson  number  Ri  veraua  particle  concentration  N  cm“3 
for  channel  2  for  the  period  06h  12.2.78  to  18h  14.2.78. 
Richardson  number  Rl  versus  particle  concentration  N  cm"3 
for  channel  5  for  the  earne  period  ea  Figure  5. 

Lou-level  radiosonde  ascents  data  at  Aughton  on  14th 
February  1978.  Dry  bulb  temperatures  only  at  OOh  (•) 
and  12h  ( - ) 

Airetream  trajectories  to  Manchester  in  the  period  21st- 
27th  February  1978. 

1  22.2.78;  2  23.2.78;  3  24.2.78 
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9  Diurnal  variation  of  particle  concentrations  N  cm 
channel  size  categories  1-5  raapectively 

10  Wind  speed  U,  relative  humidity  end  temperature  T 
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versus  pertlclr  concentration  N  cm  for  channel  1. 


Figure  11  Ulnd  epead  U,  ralatlva  humidity  and  tamparatura  T 
varaua  particle  concentration  N  cm'3  for  channel  5. 
Figure  12  Synoptic  altuation  at  12h  17th  September  1977. 

Figure  13  Airmaaa  trajectories  to  Great  Oun  Fall  for  the  period 
16th-19th  September  1977 
1  16,9.77?  2  18.9.77?  3  19.9.77 

Figure  1<*  Particle  concentration  N  cm*3  varaua  time  for  Great 
Oun  fell  from  16-19th  September  1977.  Channel  sizes 
Indicated  are  those  for  Figure  2. 

Figure  15  Variation  of  the  Inversion  height  H  with  time  from 
16th-T9th  September  1977. 

Height  A  la  the  level  of  Great  Oun  Fell  aummlt. 

Figure  16  Low-level  radiosonde  ascent  data  for  Ghanwell  on  19th 
September  1977.  Dry  bulb  temperatures  at  OOh  (— )  and 
I2h  ( - ).  The  summit  level  of  Great  Dun  Fell  la 


Indicated  (A). 
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